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ABSTRACT 

We present new measurements of the temperature-density (T — p) relation for neutral hydrogen 
in the 2.0 < z < 2.8 intergalactic medium (IGM) using a sample of ^6000 individual H I absorbers 
fitted with Voigt profiles constrained in all cases by multiple Lyman series transitions. We find model 
independent evidence for a positive correlation between the column density of H I (-/Vjji) and the 
minimum observed velocity width of absorbers (6 m in)- With minimal interpretation, this implies that 
the T — p relation in the IGM is not "inverted" , contrary to many recent studies. Fitting & m j n as a 
function of JVhi results in line width - column density dependence of the form 6 m i n = bo (VHi/A?Hi,o) r 1 
with a minimum line width at mean density (p/p = 1, Vhi.o = 10 136 cm~ 2 ) of bo — 17.9 ±0.2 km s _1 
and a power-law index of (r — 1) = 0.15 ± 0.02. Using analytic arguments, these measurements imply 
an "equation of state" for the IGM at (z) = 2.4 of the form T = To (p/p) 1 " with a temperature at 
mean density of T = [1.94 ± 0.05] x 10 4 K and a power-law index (7 - 1) = 0.46 ± 0.05. 
Subject headings: intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 

The "equation of state" of the low-density intergalac- 
tic medium is believed to be controlled by two princi- 
ple processes: photo- heating and adiabatic cooling. The 
cooling is most directly tied to the expansion of the Uni- 
verse, while the heating is expected to be a complicated 
mixture of relic effects from the reionization of hydrogen 
and helium plus the c urrent heating, predominantly from 



the U V background (Hui & Gnedin 1997 Schaye et al 
19991. This naturally imposes a relationship between 
the temperature and density of intergalactic gas. Denser 
gas is expected to trace larger over-densities, for which 
adiabatic cooling is suppressed because such regions are 
more bound against the expansion of the Universe. At 
the same time, denser gas has a much faster recombina- 
tion time scale and thus presents a larger cross-section 
for photo-ionization which leads denser gas to experience 
greater heating. 

Recently, this simple picture has been called into ques- 
tion on the basis of several statistical studies of the 



Lyman-a forest ( 


Becker et al. 2007 Bolton et al.| 


2008 


Viel et all 120091 ILidz et al.||201U| |Calura et al| 


2012 


Garziili et al. 2012|). These studies require comparison 



physical implications of the detailed shape of these flux 
distributions. Further, the shape of the distribution may 
be sensitive to the full thermal history of the gas, not just 
the present temperature-density (T — p) relation (Peeples 
et~aLp)T0" ). 
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Becker et al. (|2007|) considered the transmitted flux 
rsTribT 



probability distribution fu nction (PDF) and compared i t 



to the numerical models of Miralda-Escude et al. ( 2000 1 
finding that none of the models could explain the obser- 
vations. They found that aT-p relation in which lower 
density regions were hotter (an inverted T — p relation) 
provided a better fit to the data. Since this study, several 
other authors have obtain ed similar results using trans- 
mitted flux PDF analysis (|Bolton et aT1|2008[ |Viel et al 



2009 



Calura et al. 20121, wavelet analysis (|Lidz et aT 



2010), or a combination of both (Garziili et al 



In contrast, early work on the J' — p relation and its 
evolution with rcdshift was performed via Voigt profile 
fitting of individual H I absorption lines. These early 
studies found a monotonic T — p relation with a positive 
power-law index (|Schaye eT al. 2000; Ric otti et al. | |2000 
Bryan fc MachaceTEpUDO} [McDonald et al.|200fj K 

In this letter, we return to the more direct test of the 
T — p relation in the IGM using Voigt profile fits to indi- 
vidual H I absorbers. This method relies on the expected 
relationship between the column density of neutral hy- 
drog en of an absorber, Nhi, and i ts local overdensity, 
p/p pave et al.||1999| |Scha yc 2001). If such a relation- 
ship exists, then one expects to observe a correlation be- 
tween JVhi and the velocity widths of individual ther- 
mally broadened absorbers, b^. Absorbers having only 
thermal broadening 4 are expected to have the smallest 
bd at a given iVm, 6 m in- Thus, by observing the behavior 
of the low-6d "edge" of the distribution of &d versus Vhi, 
it is possible to infer the relationship between T and -/Vhi 
and through theory, the IGM T — p relationship. 

Here we present analysis of individual absorbers fit- 
ted in 15 high-resolution, high-S/N spectra of luminous 
QSOs at 2.5 < z < 2.9 from the Keck Baryonic Structure 

4 Other sources of broadening of absorbers include bulk motions 
of the gas (generally parametrized as a turbulent component of 
fed) and differential Hubble flow which broadens the absorption 
lines originating from the most diffuse and physically extended 
absorbers. 
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Table 1 

KBSS Absorption Line Sample 



Name 


z QSO a 




range 


S/N b 


S/N b 






Lya 


Ly/3 


Q0100+130 (PHL957) 


2.721 


2.0617- 


- 2.6838 


77 


50 


HS0105+1619 


2.652 


2.1561- 


- 2.6153 


127 


89 


Q0142-09 (UM673a) 


2.743 


2.0260- 


- 2.7060 


71 


15 


Q0207-003 (UM402) 


2.872 


2.1532- 


- 2.8339 


82 


55 


Q0449-1645 


2.684 


2.0792- 


- 2.6470 


73 


11 


Q0821+3107 


2.616 


2.1650- 


- 2.5794 


50 


33 


Q1009+29 (CSO 38) 


2.652 


2.1132- 


- 2.6031 


99 


58 


SBS1217+499 


2.704 


2.0273- 


- 2.6669 


68 


38 


HS1442+2931 


2.660 


2.0798- 


- 2.6237 


99 


47 


HS1549+1919 


2.843 


2.0926- 


- 2.8048 


173 


74 


HS1603+3820 


2.551 


2.1087- 


- 2.5066 


108 


58 


Q1623+268 (KP77) 


2.535 


2.0544- 


- 2.4999 


48 


28 


HS1700+64 


2.751 


2.0668- 


- 2.7138 


98 


42 


Q2206-199 


2.573 


2.0133- 


- 2.5373 


88 


46 


Q2343+125 


2.573 


2.0884- 


- 2.5373 


71 


45 



b The average signal to noise ratio per pixel of the QSO spectrum 
in the wavelength range pertaining to Lya and Ly/3 absorption. 



Survey (KBSS; [Rudie et al.|[2012| . We discuss the data 
set used in this study as well as the line-fitting proce- 
dure in Section [2] In Section [3j the fit to the minimum 
bd (frmin) as a function of A/hi is presented, followed by 
Section|4]in which the physical implications of the results 
in the context of the T — p relation are discussed. The 
results and their implications are summarized in Section 

Throughout this paper we assume a A-CDM cosmology 
with H = 70 km s _1 Mpc" " 



fir, 



0.3, and fi A = 0.7. 



2. DATA AND ANALYSIS 

The data analyzed in this paper are taken from the 
Keck Baryonic Structure Survey (KBSS). The KBSS is a 
large spectroscopic survey designed to study the gaseous 
surroundings of high-z star-forming galaxies by combin- 
ing redshift surveys for galaxies at 2 < z < 3 in the fields 
of bright QSOs with absorption line data from those same 
QSOs. The data include high-resolution, high signal-to- 
noise echelle spectra of 15 luminous (my ~ 15.5 — 17) 
QSOs located near the centers of the KBSS fields. The 
KBSS will be described fully by Steidel et al. (in prep); 
however, the data of relevance to this paper as well as 
the analy sis of the absorption -line sample have been pre- 
sented in |Rudie et al.| ( |2012[ ) . Here we give a brief sum- 
mary. 

The QSO spectra were obtained wi th the High Reso 



lution Echelle Spectrometer (HIRES; Vogt et al.||1994 | 
on the Keck I telescope. The HIRES spectra have 
R ~ 45,000 (FWHM- 7 km s" 1 ) and S/N - 50 - 200 
per pixel. They cover the wavelength range 3100 - 6000 
A with no spectral gaps, allowing for the observation of 
Ly/3 A1025.7 down to at least z = 2.2 in all 15 of our 
sightlines. The additional constraints provided by Ly/3 
(and in many cases, additional Lyman series transitions) 
allow for highly accurate measurements of column den- 
sities and (most importantly for this paper) line widths 
for the H I absorbers in our sample. 

The reduction of the QSO spectra and a detailed de- 
sc ription of the proce ss of fitting the forest are discussed 



inll 
full 



Rudie et al. (2012). The final Voigt profiles fits to the 



pleted using VPFIT 5 written by R.F. Carswell and J.K. 
Webb. The redshift range included in the fit is given in 
Table [T] In this paper, we use the pathlength weighted 
mean re dshift of the sample (z) = 2.37 as the fiducial 
redshift flBahcall fe Peebles||1969 |). 6 

The final H I absorber catalog includes 5758 absorbers 
with 12.0 < log(A r m/cm- 2 ) < 17.2 and 2.02 < z < 2.84 
over a total redshift path length of Az = 8.27. This H I 
sample is the largest ever compiled at these redshifts and 
is more than an order of magnitude larger than previous 
samples that included constraints from higher-order Ly- 
man lines. 

3. THE TEMPERATURE-DENSITY RELATION IN THE IGM 

The equation of state of the IGM is expected to have 
the form ^ 

T = T °{ff (1) 



whe re To is the tempera ture at the mean mass density 
(p) ( |Hui fc GnedmfT997| ). 

Under the assumption of a relatively uniform ionizing 
radiation field, a po wer-law relatio nship between A/hi and 
p is also expected. Schaye (2001) presented a model for 
the low-density IGlVL in which 'clouds' are in local hy- 
drostatic equilibrium and therefore typically have sizes 
comparabl e to the local J eans length. Employing this as- 
sumption, Schaye (2001) derived a relationship between 
AThi and local overdensity, p/p. Using updated cosmol- 
ogy and evaluating at the path-length weighted mean 
redshift of the sample: 



Pb/Pb 



10 13.6 



2/3 



r 12 

0.5 



2/3 



l + z 
3.4 



(2) 



where T4 is the gas temperature in units of 10 4 K and Ti2 
is the hydrogen photoionization rate in units of 10~ 12 s -1 
with t he no rmalization suggested by |Faucher^Gi gucre 



et al. (2008J. Assuming this scaling 7 , absorbers with 
log(A r Hi/cm -2 ) « 13.6 are expected to trace gas at the 
mean density of the universe at z — 2.4. 

Thermally broadened absorbers are also expected to 
follow a power law relation between 6d and temperature: 
bd oc T 1 / 2 . Combining with the expected T — p relation 
and the conversion between /Vhi and p, the relationship 
between 6 m i n and A/hi would be a power law of the form: 



= b 



A/i 



HI.O 



r-i 



(3) 



where 60 is the minimum line width of absorbers with 
Am = Ahi.o- With this formalism, (7 — 1) is proportional 
to (r — 1) (see eg. 



More explicitly, 



Schaye et al.||1999|) 



or pure thermal broadening: 
(2fc B T/m p ) 1 / 2 



(4) 



full Lya and Ly/3 forests of the 15 KBSS QSOs were com- 



5 http://www.ast.cam.ac.uk/~rfc/vpfit.html; © 2007 R.F. Car- 
swell, J.K. Webb 

6 The co moving pathlength, dX, was defined by Bahcall & Pec- 
|bles| | |1969| such that absorbers with constant physical size and 
comoving number density would have a contant number per dX. 

7 We note here that with T4 = 2 as we find in this paper, the TVjji 
corresponding to p is log(A r fji/cm — 2 ) = 13.7, a change of 0.1 dex. 
However, the uncertainty in other parameters (e.g. Ti2) is large 
enough, that we do not consider this small effect in this paper. 
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Figure 1. The doppler widths of absorbers (bd) versus their 
column density (Nm) for all the absorbers in the HI sample 
with relative errors in Nm and fed less than 50%. Point- 
density contours are over plotted to guide the eye in the bot- 
tom panel. Note that the minimum value of fed at each Nm 
increases with increasing Nm suggesting a normal T — p re- 
lation. 



where kg is the Boltzmann constant and m p is the mass 
of the proton. This suggests: 



log(T) = C + 21og 



km s 1 



(5) 



where 



relationship to be roughly: 8 

7 -lw3(T-l) 



(7) 



Rearranging the above equations, we expect the con- 
version between the index of &min(-^Hi) and the T — p 



In the sections that follow, we fit to the trend of the 
femin as a function of iVm for the data sample. We then 
use equations [5] and [7] to estimate the T — p relation in 
the (z) = 2.4 IGM. 

3.1. A "normal" T ~ p relation 

The measured values of fed and Nm considered in this 
letter are presented in Figure [I] Examination of Fig- 
ure [T] already suggests the main result of this paper; fo- 
cusing on the low-fed edge of the distribution, one can 
easily see that the lower envelope increases monotoni- 
cally with Nm. Such behavior is expected for a normal 
temperature-density relationship in which denser regions 
are hotter. 

3.2. Outlier Rejection 

Before fitting to the "ridge-line" of fe,„i n versus Nm, we 
exclude those absorbers that have large erro rs in their 
measured parameters. Schaye et al. (1999) suggest a 
simple algorithm for the outlier rejection in which ab- 
sorbers with relative error in fed or Nm larger than 50% 
are excluded. 9 This method is expected to primarily re- 
ject lines which originate in blends and are thus unlikely 
to have accurate line-width measurements. Additionally, 
we exclude those absorbers with &d= 8 km s~ x and &d= 
100 km s _1 . These line widths correspond to the allowed 
line-width limits input to VPFIT and as such are artifi- 
cial. The absorbers which remain after those exclusions 
are shown in Figure [Tl 

Considering Figure |ij we note a small set of points at 
very low fed which seem to lie significantly below the main 
locus of points. These absorbers are likely to be uniden- 
tified metal line contaminants, and thus we consider a 
er-rejection algorithm to exclude them. 

The cr-rejection algorithm is applied as follows. Only 
absorbers with fed < 40 km s _1 are considered. The data 
are sorted by their Nm- The absorbers are placed into 
./Vhi bins of width 0.25 dex and the mean and standard 
deviation of fed are computed for the bin. Absorbers with 
fed 2(7 or larger from the mean are flagged and excluded. 
This process is iterated; the set of absorbers excluded via 
this method with &d lower than the mean are surrounded 
by (blue) boxes in Figure [2] and are excluded from the 
fitting procedure. We have checked the majority of a- 
rejected absorbers by hand and find that they are either 
previously unidentified metal lines or parts of blended 
saturated H I systems whose individual properties are 
poorly constrained. 

In the following section, we proceed with the measure- 
ment of the femin - Nm relation. We show that the ap- 
plication of the cr-rejection algorithm does not produce 
significant changes to the results. In the conclusions of 
the paper, however, we prefer the combined error and cr- 
rejection method as it excludes the majority of the poorly 
measured absorbers as well as those which lie far below 
the & m i n - Nm trend. 

8 In the following equation, we have neglected the minor depen- 
dence of Pb/J>b on ^ from equation [2] 

9 The relative error considered is computed from the formal error 
bars output by VPFIT for both b d and Nm- 
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Figure 2. Left panel: bd versus Nm for all absorbers in the KBSS HI sample with 12.5 < log(A r Hi/cm -2 ) < 14.5 and relative 
errors in 6d and Nhi less than 50%(black points). The curves show the resulting best power-law fits deter mined by applying th e 
bootstrap resampling method to the sample followed by the iterative power-law fit method described by Schaye et al. (19991. 
Top panels: The blue and red lines pertain to the data set with the cr-rejected absorbers (blue boxes) excluded or included 
respectively. Top Center and Right: The probability distribution for the fitted parameters using the 50% relative-error rejection 
algorithm (red) and both the 50% relative-error rejection algorithm and the 2a rejection algorithm (blue). Bottom panels: The 
orange and green curves pertain to the z > 2.6 and z < 2.6 subsamples of the cr-rejection data set. Here JVhi.o varies between 
the subsamples with log(A r Hi,o/cm -2 ) = 13.7 for z < 2.6 and log(A r Hi,o/cm~ 2 ) = 13.4 for z > 2.6 due to the redshift dependence 
in equation [2] Note that the bmm-A^Hi relation is roughly constant, but the expected AThi,o evolution would require evolution in 
feo and To- 



3.3. Fitting b m - m versus Nm 

In this section we measure the 6 m i n - Ahi relation 
using the two data sets resulting from the outlier re- 
jection methods described above. To fit the trend be- 
tween b m \ n and Nm, we follow the i terativ e power-law 

A power- 



(1999). 



method proposed by Schaye et al. 
law relationship of the form shown in equation [3] with 
log(A r Hi,o/cm _2 )= 13.6 is fit to the data. Data points 
which have 6d larger than one mean absolute deviation 
above the fit are discarded. The power law is refit to the 
remaining absorbers and the rejection and refitting are 
repeated until the power law converges. Then, absorbers 
more than one mean absolute deviation below the fit are 
removed once, and the power law is refit. This fit is taken 
as the final form of the minimum 6 m ; n - Nm relation. Er- 
rors in the parameters of the fit are derived by applying 
the bootstrap resampling method to the data sample and 
following the above outlined procedure. 

Figure [2] shows the results of the iterative power-law 
fit for both data sets formed via the outlier rejection al- 



gorithms. The red curves correspond to the fit including 
the CT-rcjectcd points, while the blue curves correspond to 
those without. The center and right panels show the dis- 
tribution of fit parameters obtained through bootstrap 
resampling. As expected, those fits which exclude the 
low-&d outliers via cr-rejection (blue histograms) result in 
tighter distributions for the fitted parameters; however, 
the best fit values are nearly indistinguishable. 

The best fit parameters bo and (T — 1) and their rela- 
tive errors, as well as their physical counterparts To and 
(7 — 1), are presented in Table [2] To and (7 — 1) are 
calculated assuming equations [5] and [7] respectively. No- 
tably, both outlier rejection algorithms produce statisti- 
cally consistent results suggesting that the measurements 
are robust to changes in the rejection algorithm; however, 
(as expected) slightly higher b are preferred by the al- 
gorithms with cr-rejection. 

Henceforth, we consider the results from the data set 
formed with both the error rejection and cr-rejection al- 
gorithms applied prior to the power-law fit to & m in(-^Hi)- 
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Table 2 

Fits to the fe m i n — Afjji and T — p relation in IGM 



Outlier Rej. z range (z) log(AT H i,o/cm- 2 ) (T - 1) (7 - l) a b T /10 4 K b 



default 2.0 - 2.8 2.4 13.6 0.156 ± 0.032 0.47 ± 0.10 17.56 ± 0.40 1.87 ± 0.08 

2cr-rej. 2.0 - 2.8 2.4 13.6 0.152 ± 0.015 0.46 ± 0.05 17.90 ± 0.21 1.94 ± 0.05 



2<x-rej. 2.0 - 2.6 2.3 13.7 0.156 ± 0.016 0.47 ± 0.05 18.50 ± 0.22 2.07 ± 0.05 

2<j-rej. 2.6- 2.8 2.7 13.4 0.171 ± 0.032 0.51 ± 0.10 17.39 ± 0.44 1.83 ± 0.09 

a (7 — 1) is calculated from the measured value of (T — 1) using equation [7] 
Tq is calculated from the measured value of bo using equation I5I 



Table 3 

Comparison with Previous measurements of the T — p relationship 



Source 
|Schaye et ah] ( |2000| | 

|Ricotti et"al~] \2000\ 



Bryan fc Machacekl mhb 
McDonald et aEl 
Z;aidarriaga et ai. 
Becker et al.J |2U0'7 




|Bolton et ah| |2008[ ) c 

|Viel et al.| ( [2009} 
|Lidz et al.j | |2010[ | 
|Becker et aT] | |2011| | 



Garzilli et al. 1 2012 1 



|Calura et al] | |2012| 



Method 

Line Fitting 
Line Fitting 

Line Fitting 
Line Fitting 

Line fitting 
Line Fitting 

Transmission Power Spectra 
Flux PDF 

Flux PDF 
Flux PDF 
Flux PDF 

Flux PDF 

Wavelet Analysis 
Wavelet Analysis 

Curvature Analysis 

Wavelet Analysis + Flux PDF 
Wavelet Analysis 
Flux PDF 

Wavelet Analysis + Flux PDF 
Wavelet Analysis 
Flux PDF 

Flux PDF 
Flux PDFS 



2.3 
2.5 

2.75 
1.90 

2.7 
2.41 
2.4 

2-6 

2.07 
2.52 
2.94 



2.6 
2.2 

2.4 

2.1 
2.1 
2.1 
2.5 
2.5 
2.5 

3 
3 



T /10 4 K 

1.17 
1.38 



2.52 
1.77 

1.65 a 
1.74 ± 0.19 

1-3 



1.9 ±0.6 

1.6 ± .6 d 
2.1 ± .7 d 

1.11 ± 0.06 c 

1.7 ±0.2 
1.6 ±0.5 

1.5 ±0.3 

1.3 ±0.4 

1.6 ±0.4 

1.4 ± 0.9 

1.93 ±0.48 
1.79 ±0.35 



(7-1) 

0.27 
0.56 



0.22 ±0.10 
0.32 ±0.30 

0.29 
0.52 ±0.14 
to 0.6 
» -0.5 

-0.33 
-0.46 
-0.56 

-0.25 ±0.21 



0.11 ± 0.11 

> -0.14 
-0.01 ±0.14 

> -0.05 

> -0.08 

> -0.31 

-0.10 ±0.21 
-0.30 ±0.12 



Comment 
Their Figure 6 

Their Table 6 



Their Table 5, z slm = 3 
Their Figure 3 



Their Figure 30 
Their Table 3 

Their Table 2 
Their Table 5 



This Work 



Line fitting 



2.37 1.94 ±0.05 



Results from data set 
0.46 ± 0.05 using both relative error 
and cr-rejection 



eludes 18 QSOs (metals removed by hand to eliminate 



a The authors suggest this is an upper limit. 
b Based on th e spectrum of HS 1946+7658 
c Measured by |Kim et al.| (2007) , sample 
d Marginalized over (7 — 1) — — 0.6 
e Assuming 7 — 1 — 0.56 

^ Find that flux PDF results in lower v alues of 7 than wavelet analysis in general. 
s Joint analysis with [Kim et al~]j2007| . 



noise in 



flux PDF.) 



4. RESULTS 

As discussed above, the mapping of 6 m i n (JVHi) to the 
equation of stat e for intergalactic ga s requires a theoret- 



ical framework. |Schaye et al.| ( 2000 ) considered simula- 
tions of the IGM and found that the minimum values of 
bd at each iVui were larger than expected given purely 
thermal broadening. They suggested that all absorbers 
in the Nm range used to measure the T — p relationship 
experienced mild Hubble broadening (also expected to 



scale as T 1 / 2 for clouds having sizes similar to the local 
Jeans scale.) In this case, estimates made using the above 
equations would mildly over-predict the temperature at 
mean density. We prefer to report physical values tied 
closely to our measurements rather than rely directly on 
a specific set of simulations; however, our measurements 
of 60 and r — 1 can be converted to physical parameters 
using a different set of assumptions. Our simplifying as- 
sumption of pure thermal line broadening for absorbers 
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near b m i n (Nm) results in Tq values slightly higher than 
some previous studies which used simulations as a refer- 
ence (see Table [3]). 

In Table |3j we compare our results with those of previ- 
ous studies. We find generally good agreement between 
our results and those of past studies which used line fit- 
ting. We further emphasize that our results differ signif- 
icantly from those obtained using the transmitted flux 
PDF and similar statistical methods, calling into ques- 
tion the interpretation of such techniques. 

4.1. He II Reionization 

One of the expected results of He II reionization, be- 
lieved to occur at z sa 3, is a significant change in the 
T—p relation du e to excess photoionization heating ( Hui 
fc Gnedin||1997 1 . These changes are not long lived, and 
the exact effects depend on the speed and patchiness with 
which reionization procee ds and the spectral ha rdness of 
the ionizing sources (see McQuinn et al. 2009[ ). While 
detailed discussion is beyond the scope of this letter, we 
briefly mention that the absorbers in the KBSS H I sam- 
ple show no strong evidence for a change in the slope 
of bmin(-^Hi) for higher-redshift absorbers, which (pre- 
sumably) are temporally closer to the He II reionization 
epoch. The bottom panels of Figure [2] show the iterative 
power law fit to absorbers with, respectively, z < 2.6 
(green) and z > 2.6 (orange) from the sample with a- 
rejection. Notably, the 6 m i n — AThi relation appears to 
be independent of redshift (left panel); however, the ex- 
pected evolution in the value of Nm at the mean density 
(equation p2|) results in differing values of AThi.o f° r the 
two redshift bins. This in turn leads to expected vari- 
ation in the values of bo and To with z. The implied 
parameters from the fit are listed in Table [2] 

5. CONCLUSIONS 

We have inferred the T — p relationship in the (z) = 2.4 
IGM using Voigt profile fitting of individual Lyman 
lines. Fitting the trend of the minimum line width 
(frmin) as a function of column density (AThi) with a 
power law, we find best fit values of 6 m i n at mean den- 
sity (which at (z) = 2.4 corresponds to AThi.o = 10 13 ' 6 
cm~ 2 ) 60 = 17.9 ± 0.2 km s _1 and a power-law index 
(r — 1) = 0.15 ± 0.02. Assuming a monotonic relation 
between AThi and p/p, these data support the conclu- 
sion that the temperature-density (T — p) relationship 
in the IGM is not inverted at (z) = 2.4 but instead 
has a power-law index (7 — 1) = 0.46 ± 0.05. Further, 
our results suggest a temperature at mean density of 
T = [1.94 ± 0.05] x 10 4 K. Within our sample spanning 
the redshift range 2.0 < z < 2.8, there is no evidence for 
significant evolution in the 6 m ; n — Nm relation or in the 
power-law index (7 — 1) of the T—p relation. 
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